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ABSTRACT 

Recent gravitational microlensing observations predict a vast population of free- 
floating giant planets that outnumbers main sequence stars almost twofold. A 
frequently-invoked mechanism for generating this population is a dynamical in- 
stability that incites planet-planet scattering and the ejection of one or more 
planets in isolated main sequence planetary systems. Here, we demonstrate that 
this process alone probably cannot represent the sole source of these galactic 
wanderers. By using straightforward quantitative arguments and N-body simu- 
lations, we argue that the observed number of exoplanets exceeds the plausible 
number of ejected planets per system from scattering. Thus, other potential 
sources of free-floaters, such as planetary stripping in stellar clusters and post- 
main-sequence ejection, must be considered. 
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1 INTRODUCTION 

One possi ble explanation f or the e xistence of free-floating 
planets dLucas fc Roche] l2000l: IZapatero Osorio et al] 
I2000L 120021 : iBihain et al.M2009l : ISumi et al.ll201lf ) is that 
they formed in protoplanetary disks around young stars, 
in systems of multiple planets. These planetary sys- 
tems subsequently underwent large-scale dynamical in- 
stabilities involving close encounters between planets 
and strong planet-planet scattering events that ejected 
some pl anets and left t he su r viving planets on perturbed 
orbit s dRasio fc Ford 19961: Weidenschilling fc Marzaril 
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tribution of orbital separations b etween adjacent two- 
planet pairs (iRavmond et al.ll2009l), and perhaps certain 
resonant systems (|Ravmond et alJ ^OOS). 

In order for planet-planet scattering to create the 
free-floating planet population, the following equation: 



N s tar 



— f giant X /unstable X 7l e jec, 



(1) 



must hold, where Nf ree /N s t ar s is the observed frequency 
of fre e-floating giant p lanets of 1.8lJ'g per main-sequence 
star |Sumi et al.ll20lj) , f giant is the fraction of stars with 
giant planets, j "unstable is the fraction of giant planet sys- 
tems that become unstable, and the mean number 
of planets that are ejected during a dynamical instability. 
The terms on the right-hand side of Eq. {TJ are all depen- 
dent on stellar mass. We dis cuss these correlations exte n- 
sively in Section 3; see also lKennedv fc Kenv on (2008). 

Exoplanet observations constrain the frac- 
tion of stars w ith gas giant p lanets to be larger 
than ~ 14% (ICumming et all |200S| ; iHoward et aH 
' 

50% 
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.20 111) and perhaps as high as 
20101 ). The maj ority of giant plan- 
beyond 1 AU l|Butler et al l 120061 : 
7]j, and their abundance increases 
distance within the observa- 



ijGould et al] 

are locate d 

lUdrv fc Santos! 12007 
strongly with orbital 
tional capabilities (~ 5 AU) of radial velocity surveys 
(|Mavor et alj|201ll ). 

Given the difficulty of me asuring eccentricities 
with radial velocity m easurements (jShen fc Turnerll2008l : 
IZakamska et al.l l201ll ). the fraction of planetary sys- 
tems that becomes unstable is modestly constrained 
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by the eccentricities of surviving planets. In addi- 
tion, there is a clear positive mass-eccentricity corre- 
latio n: more massive exoplanets have higher eccentrici- 
ties jRibas fe Miralda-Escudell2007l ; lFord fc Rasioll200Sl ; 
IWright et al.l 20091 ). The simplest way to reproduce the 
observed distributions is if a large fraction of systems 

- at least 50% but more probably 75% or more - be- 
come unstable, and if the giant planets' masses within 
system s with high-mas s (M > Mj) planets are roughly 
equal (Ra ymond et al.l I2010T ). The typical number of 
planets ejected per unstable system, n e j ec , must be an 
increasing function of the number of planets that form 
in a given system. However, n e j ec has been addressed 
only tangentially in previous studies, and we quantify 
this value in a consistent manner here. 

Assuming observationally motivated constraints - 

N free/ N stars = 1-8, fgiant = 0.2, and funstable = 0.75 

- each instability must eject 12 Jupiter-mass planets 
(i.e., n e jee = 12). For the full range of plausible con- 
straints - N free /Nstars = 1 - 3.5, fgiant = 0.14 - 0.5, 

and /unstable = 0.5 — 1 - the range of values for n e jee 
is between 2 and 50. At first glance, these values appear 
implausibly large. And as we will show, n e jec is indeed so 
large as to be inconsistent with observational constraints, 
which means that the planet-planet scattering model, at 
least in its simplest form, cannot explain the free-floating 
planet population. 

Section 2 is dedicated to placing theoretical con- 
straints on n e jec using N-body simulations of planet- 
planet scattering. We chose initial conditions that we 
expect to be the most efficient at ejecting planets and 
tested systems initialized with up to 50 planets. Section 
3 incorporates this result into our main argument and dis- 
cusses additional considerations, such as the contribution 
from previous generations of main sequence stars which 
are now stellar remnants, and other potential sources of 
free- floaters. 



2 SCATTERING SIMULATIONS 

The dynamical evolution of unstable multi-planet sys- 
tems leads to planet-planet collisions, planet-star colli- 
sions, and/or the hyperbolic ejection of one or more plan- 
ets. Previous work has s hown that ejections comprise the 
most frequent outcome fWeidenschilling fc Marzarill 19961; 



Papaloizou fc Terquemll200ll; iMarzari fc Weidensc hilline 
20021; lAdams fc Laughlinl 120031; lJuric fc Tremainel |200S| ; 
Veras et all 120091 : iRavmond et al.1 I2010T I. In an exten- 
sive investigation of pla net-planet scatteri n g in s ystems 
of 10 and 50 planets, lJuric fc Tremainel (|2008T ) found 
that 50% — 60% of all planets were ejected from unstable 
systems with little depe ndence on the in i tial p lanetary 
distribution. In addition, IRavmond et alj (|2010l ) showed 
that the number of ejected planets is larger in systems 
that start with unequal-mass planets and in systems with 
higher-mass planets. 

In order to self-consistently assess n e j cc as a func- 
tion of the number of planets which have formed 
in a planetary system, N p , we performed a suite of 
planet-planet scattering simulations. Although the phase 
space of initial conditions is extensive, the results of 



the investigations refe renced above (particularly from 
lJuric fc Tremainell2008l ) demonstrate that n c j cc is largely 
insensitive to the initial distribution of planetary eccen- 
tricity or inclination. Therefore, we consider initially cir- 
cular and nearly-coplanar (i chosen randomly from 0° — 
1°) sets of N p = 3 — 50 planets. We tested two different 
planetary mass distributions: in the equal simulations 
each planet was 1 Jupiter-mass, and in the random simu- 
lations, the logarithm of the mass of each planet was cho- 
sen randomly such that the planet masses took values be- 
tween 1 Saturn-mass and 10 Jupiter-masses. We assumed 
that giant planets tend to form at a few AU - close to the 
ice line - in marginally-unstable configurations. Thus, the 
innermost planet was placed at 3 AU and for N p — 3, 4, 5, 
6, 8, and 10, additional planets were spread out radially 
with a separation of K = 4 mutual Hill radii between ad- 
jacent planets (as K determines t he instability timescale ; 
IMarzari fc Weidenschillindl2002l ; IChatteriee et~al1l2008l ). 
For N p — 20 and 50, K was decreased to keep the out- 
ermost planet at < 200 AU (for N p = 20 and 50, K 
= 2.56 and 0.99 for the equal simulations, and K — 
2.00 and 0.70 for the random simulations). Each system 
was integrated for 10 Myr using the Bulirsc h-Stoer in- 
tegra tor in the Mercury integration package (|Chambersl 
1999). Although slow, this integrator accurately models 
close encounters. The radius of each planet was taken to 
be Jupiter's current radius, and collisions were treated as 
inelastic mergers. A planet was considered to be ejected if 
its orbital distance exceeded 10 5 AU, which represents a 
typical distance at which galactic tides can cause escap e 
over a typical main sequence lifetime (|Tremainel 1 19931 ). 
We ensured that each simulation conserved energy and 
angular momentum to one part in 10 4 or b etter to avoid 
numerical artifacts (I Barnes fc Quinnll200l ). 

We obtained a sufficiently large statistical measure 
of n e jec by running ensembles of systems for each value of 
Np. We integrated 100 sets of initial conditions for each 
of N p — 3, 4, 5, 6, 8, and 10, and assigned a random 
value to each planet's orbital angles in each instance. We 
did the same for N p = 20 and 50, but instead ran 40 and 
10 simulations, respectively, for each due to the increased 
computational cost. Each simulation was run for 10 7 yr. 
Our results are displayed in Fig. 1. 

The figure demonstrates that between 20% and 70% 
of planets are ejected in each set of simulations. The ver- 
tical bars attached to the mean escape percentages for 
each bin represent ±1 standard deviation values from the 
mean. For the equal simulations (left panel) this mean 
ejection fraction is roughly constant with N p , and varies 
by at most 10%. In contrast, the ejection rate in the 
random simulations (right panel) increases monotonically 
with Np. The mean ejection percentage for N p — 10 is 
50%-60%, cor roborating the res u lts fr om the 10-planet 
simulations in ljuric fc Tremainel l|2008l ). The red dashed 
lines represent the fraction of the initial planetary mass 
ejected for each value of N p (slightly offset to the right 
of the blue curves for clarity). This value is less than 
the percentage of planets ejected in all cases, demon- 
strating that the most massive planets tend to scatter 
out the smaller o nes, corroborating previous results (e.g. 
iFord et al.ll2003l ). 

The effect of varying the initial system configura- 
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Figure 1. Ejection fractions for planetary systems with different initial numbers of giant planets (N p ). The left panel shows the 
results for the equal simulations, in which all planets have the same mass (1 Jupiter- mass). The right panel shows the results 
of the random simulations, where the planetary masses were selected randomly and logarithmically between 1 Saturn-mass and 
10 Jupiter-masses. Vertical bars about the mean represent one standard deviation. The solid blue curves indicate the fraction 
of planets that is ejected, and the red dashed curve (which is slightly offset to the right from the blue curve in the right panel 
for clarity) represents the fraction of initial planetary mass that is lost. The four other segments in the right panel all represent 
different initial distributions for N p = 5 (green: each planet has eo = 0.5; pink: each planet's initial eccentricity is assigned a 
random value; orange: M* = 3Mq and a,i„ nermoat = 30 AU; purple: = (1/3)Mq and dinnermost = 0.3 AU). These plots 
provide relations between N p and n ejec , the number of planets ejected; values from the unequal mass plot are used as coefficients 
in Eq. J2}. 



tions does not have a drastic effect on t he ou tcome, as 
demonstrated in, e.g. Jjuric fc Tremainel (|2008T l. In order 
to sample this effect for our setup, we have performed 
4 additional sets of simulations, all with N p = 5, and 
have plotted the results on Fig. 1 to the immediate left 
of the N p = 5 blue curves. The green curves demonstrate 
the ejection percentage when each planet is initialized 
with e = 0.5 and the pink curves for a randomized dis- 
tribution of initial eccentricities. The orange and purple 
curves assume respectively different stellar masses (3M© 
and O.3A/0) and different innermost semimajor axes val- 
ues (0.3 AU and 30 AU) in order to account for the dif- 
ferent approximate ice line boundaries. In all these cases, 
the fraction of mass ejected is less than the fraction of 
planets ejected, as in the nominal N p = 5 cases, and the 
mean ejection fraction of planets varies by about 20%. 



3 DISCUSSION 

We have related Nf ree / N s tars to n e jec in Section 1, and 
have related n e jec to N p in Section 2 and Fig. 1. Now, 
we can link the two relations, and do so through a more 
precise formulation of Eq. (TjTJ): 



Nfr 



f (N P ) „ ;( n p) v „( N P ) 

able A n ejee 



/ j J giant J unstab 
iV p =3 

f(3) f(3) .ntv.A.fW tW 



— f) 1 7 ■ S ■ f { ' f { ' 4- "V7 ■ 4 • f y > f 1 ' 
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+... (2) 

where the coefficients are taken from the mean ejection 



fractions of the random simulations in Fig. 1. and which 
is subject to the constraint 



giant 



< 1. 



(3) 



Now we can simply set N free/ N a tars = 1.8 and take a 
closer look at Eqs. ©-Q. 

In order to satisfy the constraint in Eq. ((3}, the vast 
majority of planetary systems must have N p ^ 5. If wc 
use the observation ally determined upper bound of 0.5 
(jGould et alJ l2010h on the right-hand side of Eq. Q, 
then the vast majority of planetary syste ms must have 
N p > 7. If we use the lower bound of 0.14 |Mavor et al.l 
l201ll ) - which applies for Sun-like stars - then the vast 
majority of planetary systems must have N p 20. These 
minimum values of N p would be even higher if we in- 
stead derived Eq. @ from the equal simulation results. 
If, motivated by the monotonicity of the coefficients in 
Eq. we assume that all systems form a fixed number 
of planets, then we can calculate a critical giant planet 
frequency fg^ a P „\ crit below which the free-floating popula- 
tion cannot be reproduced. We plot this quantity in Fig. 



Here, coefficients 



[2] for two different values of f 
for values of N p that were not sampled in our simula- 
tions were conservatively fixed at the coefficient of the 
next highest value of N p that was sampled. For exam- 
ple, the coefficients of fg*° nt and f giant are 0-65 and 0.70; 
(21) 

hence, the coefficient of f giant is taken to be 0.70. 

Now that we have determined the relation between 
N p and fgf^t ■ > we can assess the constraints on form- 
ing N p giant planets in a single system. A simple first 
consideration is the outermost planet's semimajor axis 
before any scattering occurs. For simplicity, assume N p 
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Figure 2. The critical fraction of stellar systems that must 
initially contain at least N p giant planets in order to pro- 
duce the free-floating planet population. The curves were ob- 
tained from Eq. assuming N free/ N s tars = 1.8^J'g; the 
solid curves with symbols assume the nominal value of 1.8, 
and the dashed curves without the symbols assume values of 
1.0 and 3.5. All three blue curves assume every system be- 
came unstable, and the three red curves assume that half of 
all systems became unstable. Current observational limits for 
systems which contain any number of detected planets are 
given by the horizontal black lines. This plot illustrates that 
an unrealistically large number of giant planets must form per 
star if planet-planet scattering represents the sole source of 
the free-floating planet population. 



1- Jupiter-mass planets all orbit a 1M* star, and that 
the planets are formed close enough to be on the verge 
of instability (K — 4). Then, with an adopted in- 
nermost semimajor axis of 3 AU, a out ermost = 3 AU 
•(l^lS)^ 1 . Therefore, for N p > 13, a outermost > 200 
AU. Core accretion cannot form planets beyond ~ 35 
AU dDodson-Robinson et al] 12009). gravitational insta- 
bility has not yet be en demons trated to produce plan- 
ets beyond 200 AU (|Bossl l2006h . and only in the most 
extreme cases may planet-disc interactions cause plan- 
ets t o migrate outward beyond 200 AU (|Crida et al.l 
2009|). If giant planets were all formed by core accre- 
tion, then the a < 35 AU restriction implies N p ^ 8, 
implying that at least 40% of systems, all with 8 giant 
planets, all must have become unstable in order to pro- 
duce the free-floating planet population by planet-planet 
scattering alone. Alternatively, if giant planets were all 
formed by gravitational instability, even in extende d 
"maximum-mass" discs (|Dodson-Robinson et al.ll2009l ). 
then not more than a few giant -planet-mass clumps form 
(|Bolevll2009l ; iBolev et alj|2010h . 

Radial velocity observations suggest that, within a 



higher-mass stars (Ijohnson et al.l 12007 


: lLovis & MavorJ 


2007; [Bowler et alj|2010t 1 Johnson et al. 


2010). However, 



low-ma ss stars outnumber hi gher-mass stars by a large 
factor l)Parravano et all 1201 ll . and references therein). If 
fgiant does indeed decrease as M* decreases, then ex- 
plaining free-floaters via planet-planet scattering proves 
to be more difficult simply because fewer total stars would 
form giant planets. Consider, for example, the observed 



lower-bound frequencies of giant planets for different stel- 
lar masses deduced from radial velo city surveys for thre e 
categories of stars: M dwarfs (3%. jBonfils et alj 120111 ). 
Solar-like st ars: (14%, iMayor et all l201ll ). and A-tvpe 



stars (20%. Ijohnson et al.l lipid ). The total integrated 
value of fgiant depends on the mass ranges adopted for 
these stellar types and the particular initial mass function 
assumed. At minimum, the relative frequency of stars in 
these three bins should differ by a factor of a few (e.g., 
in the ratio 4:2:1), which yields an integrated fgiant of 
< 8.5%. If, however, there exists at least one order of 
magnitude more M dwarfs than Solar-type stars, then the 
integrated fgiant must be < 4.5%. If this value is adopted 
as the true value of fgiant , then Fig. 2 demonstrates that 
giant planet systems must be extremely crowded, with at 
minimum 25 giant planets per system. This value is al- 
most certainly unrealistic given the arguments presented 
above and current observational constraints. However, we 
note that microlensing surveys suggest that giant plan- 
ets on more distant orb its may actually b e very common 
around low-mass stars (|Gould et al.ll2010T i. in which case 
a typical giant planet system need only contain 4-10 plan- 
ets (Fig. 2). 

To what degree do stars that have already evolved off 
the main sequence contribute to the free-floating planet 
population? In order to estimate this contribution, we 
need only count the current population of white dwarfs, 
as stars that underwent supernovae or/and formed a 
black hole comprise a negligib le fraction (< 1% ) of the 
total stellar population (jParravano et all l201lh ■ One of 
the most complete and least-biased samples of white 
dwarfs represents the local population (within 20 pc of 
the Sun), which is esti mated to have a spac e density of 
4.8 ± 0.5 x 10" 3 pc" 3 (|Holberg et all 120081 ). The space 
density of stars in the Ga lactic Disc is ~ 6 x 10 _1 pc -3 
l|Binnev fc Tremainell2008l . Pg. 3), implying that the frac- 
tion of the free-floating population which has arisen from 
dynamical instability soon after formation in already- 
dead stars is on the order of 1%. This value is not great 
enough to change the results here unless N p is orders of 
magnitude higher for high-mass stars than for the lower- 
mass stars that dominate the current galactic population. 

If the planet-planet scattering model can reproduce 
the free-floating planet population in a realistic frame- 
work, then the following observational predictions must 
hold. First, the frequency of giant planets around low- 
mass stars must be high (~ 50%), in agreement with 



ment with radial velocity surveys ( 


Johnson et all 120071; 


Lovis & Mavorl 20071: Bowler et al.l 


20ld; Johnson et al.l 


20ld; Bonfils et al. 2011). This scenario could be ex- 



plained if giant planets around M dwarfs are simply too 
distant to be sampled by current radial velocity surveys, 
further implying that giant planets are found at larger 
orbital separations around low-mass stars than Sun-like 
stars. Second, systems containing a large number of gi- 
ant planets and extending to large orbital radii must be 
abundant, at least at early times before they become un- 
stable. Such systems may be detectable by direct imag ing 
of y oung stars (e.g. HR 8 799, lMarois et al.ll2008l .l 20ld . see 
also lDelorme et afll201ll ). Third, the free-flo ating planet 
population of N free/ 'N stars = 1.8 reported bv lSumi et alj 
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(|201 ll ) must be somewhat overestimated. Fig. 2 shows 
that observations are far easier to reproduce if this value 
is closer to 1. Finally, the vast majority of giant planet 
systems must go unstable. This statement, in turn, im- 
plies that giant planets on highly-eccentric orbits should 
continue to be common and that both terrestrial planets 
and debris disks should be anti-correlat e d with eccentric 
giant planets (|Veras fc Armitagel 120061 ; iRavmond et al.l 

Eml). 

Other potential sources of free-floating planets exist. 
One such source is dynamical ejection from multiple-star 
systems; on ly about two-thirds of all stars are single stars 
(lLadall2006l ) so multiple-star systems may provide a com- 
parable contribution to free-floaters. Another potential 
source arises from external forces such as passing stars, 
galactic tides, or - most relevant for close (< 100 AU) 
planets - perturbations while stars are still in their birth 
clusters. Indeed, studie s of fl ybys on planetary systems 
(jMalmberg et all I2007I . | 201 ll ) indicate these perturba- 
tions may eject planets. Recent cluster simulations have 
shown that planets with orbital radii of 5-30 AU are likely 
to disrupted by close stellar pas sages within their birt h 
clusters at a rate of up to ~ 10% (|Parker fc Quanzll201ll ). 
The cluster disruption rate is much hig her than the the 
long- term effects from passin g stars (e. g. Weinberg et all 



Il987l ) or galactic tides (e.g. lTremainelll993h . which are 
thought to cause disruption for semimajor axes of ~ 
10 s AU. Free-floating Jupiter-mass objects can also be 
forme d by collisions between high-mass protoplanetary 
discs |Lin et al.ll 19981) . but extremely dense stellar condi- 
tions are needed to cause such collisions. Other potential 
sources of free-floaters may arise from dynamical ejection 
pro mpted by mass loss during post-main sequence evolu- 
tion. jVeraiTeFlil] (|201ll ) show that planets at several hun- 
dred AU can be dynamically ejected due to mass loss from 
stars with progenitor masses greater than ~ 2Mq; these 
planets may reach such wide distances through scatter- 
ing over the main sequence lifetime of the system. The 
interaction of planetary systems with localized Galactic 
phenomena such as the tidal streams, passage into and 
out of spiral arms from radial stellar migration, Lindblad 
resonances with the bar, and passing molecular clouds 
represent other, largely unexplored, potential causes for 
ejection. Also, we cannot rule out the possibility that 
Jupiter-mass objects could simply represent the low-mass 
tail of the stellar initial mass function, and that free- 
floaters could help constrain the presence or e xtent of a 
planet-star gap in the initial mass function (e.g. IChabrierl 
120031 ; iParravano et al.ll201ll ). 
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